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Abstract. The high frequency vibration analysis of a composite laminate plate structure subjected to 
impact loads was investigated by using method of energy finite element analysis (EFEA). The time 
and space averaged energy density was used as the primary variable to form the governing differential 
equations. The multilayer laminate plate is simplified to be equivalent isotropic plate using the 
average concept, such as the average damping loss factor and the average group speed. The global 
system of EFEA equations can be solved numerically and the energy density distribution within the 
whole system can then be obtained. The EFEA numerical results for composite laminate plate 
structure are validated through comparison with those of very dense conventional finite element 
analysis (FEA). 


Introduction 


Fiber-reinforced composite materials have been widely used in high-performance structures that 
need to be lightweight, yet strong enough to take harsh loading conditions, such as car, ship, and 
spacecraft [1]. The innovative designs and the increasingly demanding performance requirements 
make it very importance to address the vibration issues in the design phase. 

To improve the performance of composite materials and utilize their advantages significantly, it is 
crucial to have an effective prediction of their response subjected to impact loads in practical 
engineering problems. The conventional FEA requires a large number of elements to capture the high 
frequency characteristics of the structures, which results in very long computational time and high 
costs. The EFEA formulation has been developed for making an accurate prediction of vibration 
response at frequencies where conventional FEA are no longer efficient [2]. 

EFEA method is a new method for simulating vibration response at high frequencies. A power 
flow finite element analysis is developed by Nefske and Sung [3]. Bouthier and Bernhard[4] derived 
the equations of time and space averaged energy density and intensity in the far-field. In other 
research, the EFEA formulation was applied to beams, membranes, plates and acoustic cavities[5]. 
While most of the research works of EFEA are related to isotropic materials, this paper aims to 
develop the EFEA formulation for computing the vibration response of the structures of composite 
laminate plates [6]. The global system of EFEA equations can be solved numerically and the energy 
density distribution within the entire system can be obtained with relatively small amount of mesh. 
Therefore, the vibration response of large scale composite laminate plate structures can be predicted 
in the design process efficiently. 


Energy Governing Differential Equations for Infinite Plates 


The EFEA formulation uses the time and space averaged energy density as the primary variable to 
form the governing differential equations. The expressions which govern the flow of energy in an 
infinite transversely vibrating plate were derived. The flux of energy across given closed surface, 


using the control volume approach, is equivalent to the rate of change of the total energy inside the 
surface which encloses the volume. 


Gel Ot = Pin - Paiss - V- 1. (1) 


where e is the energy density, Pj, is the input power density, Piss is the power density dissipated, V is 


the del operator, Tis energy intensity. Eq. 1 is an energy balance relationship for all elastic media and 
is valid for transient or steady analysis. 


On a time averaged basis, the potential energy and the kinetic energy density are approximately 
equal in the far-field. The energy governing equation in the infinite plate is 


- Ge / na) (0E Ox’ + FEl Oy’) + Noë = Pin. (2) 


where c, is the group speed, 7 is the damping loss factor, œ is the circular frequency. The key 
assumptions used in the derivation including harmonic excitation, steady state analysis, far-field 
approximation and symmetry about the point of excitation. A finite element formulation is developed 
to solve the Eq. 2 numerically. 


[Ke] {e} = [Oc] + [F4]. (3) 


where [Ke] is the system matrix for each element, {e} is the vector of nodal value for the time and 
space averaged energy density, [Q.] is the power flow across the element boundary, [F,] is the 
excitation vector, which represents the input energy at each node. 


Equivalent Isotropic Material for Composite Laminate Plates 


Averaged Group Speed. Group speed of transverse wave propagation for isotropic materials can be 
expressed as 


cg =2(Daw"lph)'. (4) 
where D = EWINA- Vv )], E is the Young’s modulus, A is the thickness of the plate, p is the density of 
material, vis the Poisson’s ratio. For the composite laminate plate, the group speed is a function of 


wave propagation angle a. Removing the dependency on angle a by taking the average c, from 0 to 27, 
the averaged group speed Cg can be expressed as 


E= fo" colada; e" day. (5) 


Averaged Damping Loss Factor. The total time-averaged energy loss associated with an arbitrary 
wave type can be derived using hysteretic damping model applied to each layer 


(Paiss) = Diz nwe), . (6) 


where ņ; and (e), are the damping loss factor and energy density of the /-th layer of the multi-layered 
composite plate respectively. Since (e) = (T)+(V) =2{V) , where (T) and (V) are the 
time-averaged total kinetic and potential energies, (7) and (V) are given by 


(T) = (@°/A)f, p u'udV =œ°u"Mu, (V) = (1/4) f, o'edV =u"Ku. (7) 


where u is the vector of the nodal degrees of freedom, T stands for Hermitian transpose. M and K are, 
respectively, the mass matrices and stiffness. The circular frequency for wave number k can be 
expressed in terms of the Rayleigh quotient, wo = u'K(ku / u'Mu. The damping loss factor 
associated with each type propagating wave can be expressed as 


n = Lin u'K?u/(u'Ku). (8) 


where K® is the stiffness matrix of the /-th layer. Eq. 8 is derived for each wave heading angle a;, and 
the angle averaged damping loss factor must be considered to account for the full range of wave 
propagation directions from 0 to 27. 


= fe" naddo / So" doi. (9) 


Simulation and Verification 


Using the averaged group speed and averaged damping loss factor, the EFEA differential equation 
for composite laminate plate can be expressed as 


- (C / Ha) (Cel Ax + PE Ay’) + Hoe = Pin. (10) 


From the averaged group speed C, the equivalent isotropic material property of composite 
laminate plate can be defined using equation D = C phl 16”, where p is density and A is the thickness 
of the plate. The EFEA capability for modeling and calculating the response of composite laminate 
plates in high frequencies is validated by comparing simulation results computed by EFEA with those 
of the very dense FEA. 


Fig. 1. Schematic diagram of a 4-layer composite laminate plate (30°/45°/-45°/-30°). 


The plate is a 2 m x 2 m square laminated plate with four layers of equal thickness 0.25e-3 m. The 
four layers are at 30°/45°/-45°/-30° orientations. The engineering constants of layers are 


E11=37.78e9[Pa], E22=10.0e9[Pa], E33=10.0e9[ Pa]. 
G12=G23=G31=4.91e9 [Pa], v12=0.25, p=1777.6 [kg/m?]. 


A harmonic force was applied at the center of the plate, and the input power was le-4 Watt. The 
element numbers of EFEA model are 100, while the element numbers of the very dense FEA model 
are 10,000. The computations were, respectively, carried out using the self-developed program and 
ANSYS module. The comparisons of the energy density distribution along x, y=L/2 in laminate plate 
at several frequencies are presented in Fig.2. The EFEA method eliminates the volatility of solutions. 
The numerical results of EFEA agree well with the average value of the very dense FEA, and the 
results exhibit better agreement at higher frequencies. The EFEA results can capture the energy 
distribution level in the structure well using a much smaller number of elements than that of the 
conventional FEA. The difference of the results calculated by EFEA and FEA is within 1.6dB. 


In order to analyze the complex structures, it is necessary to calculate the energy transmitted from 
the excitation location to the other components. The energy transmission coefficients for arbitrary 
joints need to be calculated and accounted for the coupling structures in the future work. 


Energy Density Distribution (f=500Hz) Energy Density Distribution (f=1000Hz) 
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Fig. 2. Distribution of energy density along x,y=L/2 in laminate plate (30°/45°/-45°/-30°) 
computed by EFEA and the very dense FEA 
(f=500Hz, 1000Hz,3000Hz,5000Hz. Reference energy densityle-12/J/m). 


Summary 


The EFEA formulation, using time and space averaged energy density as the primary variable, is 
developed for composite laminate plate structures in this paper. The composite laminate plates were 
assumed to be equivalent isotropic materials using the average concept. The EFEA method provides a 
practical approach to evaluate and predict the structural response at high frequencies, and goes 
beyond the limits of conventional FEA methods because of the computational cost or influences of 
uncertainty factors. The vibration response of the large scale composite laminate plate structures can 
be predicted in the design process and be designed to meet the increasingly demanding performance 


requirements. 
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